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ABSTRACT 

This report presents an algorithm and TI-59 
programs for position determination with Loran-C chains. 
Operational data cards are prepared in advance for Loran-C 
triplets. Position determination can be performed by 
using a single program card and an appropriate operational 
data card. 
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A. Introduction 



The Loran system is a radio aid to navigation which 
utilizes the principle of hyperbolic fixing. The locus of 
points for which the difference in arrival time of synchron- 
ized signals from a pair of transmitters is constant det- 
ermines a hyperbolic line of positions (LOP) . The 
intersection of two hyperbolic lines of position from two 
pairs of transmitters (stations) determines a position 
estimate (a fix) . That two lines of position and hence 
two pairs of stations are required for a fix does not 
necessarily mean that four separated stations are required, 
since one station of one pair may be colocated with one 
station of the other pair. In this case, a Loran triplet 
is formed (Figure 1) . Triplets may be joined "end-to-end" 
by station colocation to form a Loran chain (Figure 2). 

Loran chains are common on both the East and West Coasts 
of the North American continent. 

"Standard Loran" or "Loran-A" operates at a fre- 
quency just below 2 MHz and is now only used in the 
Pacific area. "Loran-C" operates at 100-kHz and is 
used both by ships and aircraft world wide. The compu- 
tational algorithm and programs described herein can be 
used for position determination with Loran-C triplets. 
Further information on the history, development and 
operation of the Loran systems may be found in References 
1 and 2 . 
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(a) Colocated Master Stations 




(b) Colocated Slave Stations 



M 
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{^) lolo'ited Master and Slave 



Figure 1. Loran Triplets. 
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Figure 2. Loran Chain of Five Loran Triplets. 
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B . Program Description 



The Position Determination Program requires opera- 
tional data cards (described below) for on-location position 
determination from Loran triplet time-difference measure- 
ments. Two programs are used to prepare an operational 
data card. Operational data cards should be prepared and 
validated prior to on-location navigational use. Then, 
although three programs are used in total, only one 
program is required for position determination, the 
remaining two programs are used to prepare operational 
data cards . 

The programs are : 

Master Data Card Program : This program is used to prepare 

master data cards. A master data card requires the follow- 
ing information for a master (M) station/slave (S) station 
pair (a list of current station parameters is in the 
appendix) : 

1. A M/S pair identification number. 

2. The quantity At which is the sum of the coding 
delay plus the one way base line time in microseconds 

3. The latitude and longitude of the master station. 

4. The latitude and longitude of the slave station. 

Some preprocessing of these data is performed before the master 
data card is generated. The data generated require only one 
side of a TI-59 magnetic card for each M/S pair, thus a 
second M/S pair may be placed on side 2 of a master data 

card. A master data card should be prepared in advance for 
each M/S pair that might be used within an area of operation. 
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Operational Data Card Program ; This program is used to 
prepare an operational data card for each Loran triplet 
that might be used within an area of operation. Each 
operational data card contains data merged from two 
master data cards that contain the M/S pair information 
for each pair of a triplet. These merged data are validity 
checked, colocation of master or slave determined and 
encoded . 

The inputs required to generate an operational data 
card for a Loran triplet using this program are those 
contained in the two master data cards that are associated 
with the triplet. It is possible to prepare and store 
operational data cards rather than master data cards. This 
may be desirable if there is no scarcity of cards and storage 
space, however the number of possible Loran triplets is 
considerably larger than the number of M/S pairs. 

Position Determination Program : This program is used for 

position determination. The required inputs for the program 
are the indicated time differences T for the two M/S 
pairs of the triplet and the data contained in the opera- 
tional data card associated with the triplet. The outputs 
of the program are the computed latitude and longitude of 
the Loran fix. Note: Every Loran determination gives 

two solutions; the proper solution and the spurious solution. 
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The spurious solution can almost always be rejected by 
inspection. However, if the stations of the Loran triplet 
are nearly aligned this may not be the case, and additional 
information will be required to determine the proper 
solution. 
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C . TI-59 Calculator Programs 



1. User Instructions 



Master Data Card Program (MDCP) 



Input Output 



Step 


Instructions 


Data/Units 


Keys 


Data/Units 


0. 


Turn the calculator on. If the cal- 
culator has been on, turn it off, 
then on to initialize it. 








1. 


Read side 1 of the MDCP program card. 








2. 


Press CLR and then read side 2 of the 
MDCP program. 




CLR 




3. 


Input a unique ID number for the Loran 
pair. See Note a. 


ID 


2nd, A* 


ID 


4. 


Input the coding delay At. 


At 


2nd,B' 


At 


5. 


Input the master station latitude. 
See Note b. 




A 


0 


6 . 


Input the master station longitude. 
See Note b. 




B 


0 


7. 


Input the slave station latitude. 


(|> 


C 


0 


8 . 


Input the slave station longitude. 


X 

s 


D 


0 


9. 


Run. When the calculations are 




E 





complete, 4 will appear in the 
display; then the display will blank. 

After the display has blanked, enter 
a card in the read/write slot to 
record the master data card for the 
Loran pair. When the recording is 
complete, 4 will reappear in the 
display. See Note c. 

Notes: a. Loran pairs are coded on the navigation maps using 

designators such as 9930X, 9930Y, 9930Z and 9930W. 

It is suggested that the ID*s for these pairs be 
coded as 9930.1, 9930.2, 9930.3 and 9930.4 respectively. 
However, any consistent scheme is acceptable. 
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b. The format for latitude and longitude data input is 

of the form: DDD.MMSSFF, where 

DDD denotes degrees 
MMM denotes minutes 
SS denotes seconds 
FF denotes hundredths of a second. 

A minus sign (-) denotes South latitude and West 
longitude. See the Appendix with regard to the 
latter sign convention. 

c. When the display has blanked after 4 has been displayed, 
the calculator is prepared to write a card side with 
the information (data) contained in bank 4. When a 
card is inserted, it will be recorded. Do not restrict 
the cards motion once the card has been grasped by the 
read/write mechanism. When recording is complete, the 
card will stop and 4 will reappear in the display. 

The card can then be removed. 
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Operational Data Card Program (ODCP) 



Step 



Instructions 



Input 

Data/Units 



Output 

Keys Data/Units 



0. Turn the calculator on. If the calculator 
has been on, turn it off, then on to 
initialize it. 

1. Read side 1 of the ODCP program card. 

2. Press CLR and then read side 2 of the ODCP 

program card. CLR 

3. Press CLR and then read a master data card 

containing the first pair of a Loran triplet. CLR 

4. Press A. A 

5. Press CLR and then read a master data card 
containing the second pair of a Loran 

triplet. CLR 

6. Press B. See Note a. B 

7. Run. When the calculations are complete, 4 

will appear in the display; then the display E 

will blank. After the display has blanked, 
enter a card in the read/write slot to 
record an operational data card. When re- 
cording is complete, 4 will reappear in the 
display. See Note b. 

Notes: a. Should the display flash, then the two master data 

cards do not correspond to a Loran triplet. Both 



the latitude and longitude of the colocated stations 
must be identical on both master data cards in order 
to successfully produce an operational data card. 

b. When the display has blanked after 4 has been displayed, 
the calculator is prepared to write a card side with 
the information (data) contained in bank 4; and, when 
a card is inserted, it will be recorded. Do not 
restrict the cards motion once the card has been 
grasped by the read/write mechanism. When recording 
is complete, the card will stop and 4 will reappear 
in the display. The card can then be removed. 
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Position Determination Program (PDP) 



Step 

0 . 

1 . 

2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 

9. 

10 . 



Input Output 

Instructions Data/Units Keys Data/Units 



Turn the calculator on. If the cal- 
culator has been on, then turn it off, 
then on to initialize it. 

Read side 1 of the PDP program card. 

Press CLR and then read side 2 of the CLR 

PDP program card. 

Press C to initialize, and then read the C 

operational data card for a Loran tri- 
plet that you are receiving. See Note 
a below. 

Input the observed time delay from the A 

first loran pair. 

Input the observed time delay from the T^ B 

second Loran pair. 

Compute the 1st solution latitude. 

Compute the 1st solution longitude. 

Compute the 2nd solution latitude. 

Compute the 2nd solution longitude. 

See Note b. 

To compute new solutions for new time 
delays, go to step 4. To compute new 
solutions for a new triplet, go to 
step 3. 

Note a. If the data card does not read (a flashing display results) , 
press CLR and repeat the step. If the card does not read 
on the 2nd try, the card should be cleaned (see the TI-59 
owners manual) or a new data card prepared. 

Note b. Once the proper solution (solution 1 or solution 2) has 
been determined for a triplet, then only step 6 (for 
solution 1) or step 8 (for solution 2) need be repeated 
after inputing new time delays in step 4 and 5. 



D Latitude 

R/S Longitude 

E Latitude 

R/S Longitude 
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2. Example Problem for Loran-C Pairs 9930X and 9930Y 



An example use of the Master Data Card Program using coding delays 
and latitude and longitude values from the Appendix: 

Input Output 

Step Instructions Data/Units Keys Data/Units 



0. Turn the calculator on. If the cal- 
culator has been on, turn it off, then 
on. This assures that it will be in 
the correct initial state. 

1. Read side 1 of the MDC program card. 

2. Press CLR, then read side 2 of the MDC 





program card. 




CLR 




3. 


Input the ID for 9930X. Next press 
2nd, then A. 


9930.1 


2nd, A* 


9930.1 


4. 


Input the coding delay At for 
9930X. Next press 2nd, then B. 


36389.66 


2nd,B* 


36389.66 


5. 


Input the master station latitude , 
then press A. 


34.034604 


A 


0 


6 . 


Input the master station longitude 
and press +/- for West, then press 

B. 


-77.544676 


B 


0 


7. 


Input the slave station latitude, 
then press C. 


46.463218 


C 


0 


8. 


Input the slave station longitude 
and press +/- for West, then press 

D. 


-53.102816 


D 


0 


9. 


Press E. After 4 has appeared in 
the display, enter an available 
card side in the read/write slot of 
the calculator. After recording is 
complete, 4 will reappear. Then re- 
move the card and label the recorded 
side "9300X MASTER.” 




E 




10. 


Input the ID for 9930Y, next press 
2nd and then A. 


9930.2 


2nd , A* 


9930.2 
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step 



Instructions 



Input 

Data/Units 



Keys 



Output 

Data/Units 



11. Input the coding delay At for 
9330Y, next press 2nd and then B. 

12. Input the master station latitude, 
then press A. 

13. Input the master station longitude 
and press +/- for West, then press 

B. 



52541.31 2nd,B* 52541.31 

34.034604 A 0 



-77.544676 B 



0 



14. Input the slave station latitude, 
then press C. 

15. Input the slave station longitude, 
and press +/- for West, then press 

D. 



41.151193 



-69.583909 



C 



D 



0 



0 



16. Press E. After 4 has appeared in E 0 

the display, enter an available 
card side (for example the second 
side of the card used in Step 9) in 
the read/write slot of the cal- 
culator. After recording is 
complete, 4 will reappear in the 
display. Then remove the card and 
label the recorded side "9930Y 
MASTER.’* The recorded data cards (or 
card) will be used in the next 
example . 
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An example use of the Operational Data Card Program 



step 

0 . 

1 . 

2 . 

3 . 

4. 

5. 

6 . 

7. 



Input 

Instructions Data/Units Keys 



Turn the calculator on. If the cal- 
culator has been on, turn it off, then 
on to initialize it. 

Read side 1 of the ODCP program card. 
Press CLR and then read side 2 of the 



ODCP Program Card. CLR 

Press CLR and then read the master 

data card for station 9330X. CLR 

Press A. A 

Press CLR and then read the master 

data card for station 9930Y. CLR 

Press B. See Note a of the ODCP 

instructions on page 9. B 

Press E. When the display has E 



blanked after 4 has appeared, insert 
an available card side in the read/ 
write slot. After recording is 
complete, 4 will reappear in the 
display. Label the recorded card 
side "9330X/9930Y OPERATIONAL DATA. 
Then, for identification, label the 
A key position on the card "9330X” 
and the B key position '*9330Y". 

This card will be used as the data 
card in the Position Determination 
Program example. The key labeling 
relates to that program. 



Output 

Data/Units 
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An example use of the Position Determination Program: 



step 



0 . 

1 . 

2 . 

3 . 

4. 

5. 

6 . 

7. 

8 . 
9 . 



Instructions 



Input 

Data/Units Keys 



You are in the North Atlantic and you 
are receiving 9930X and 9930Y and you 
wish to obtain a fix. 

Turn the calculator on. If the cal- 
culator has been on, then turn it 
off, then on to initialize it. 

Read side 1 of the PDP program card. 

Press CLR and then read side 2 of the CLR 

PDP program card. 

Press C to initialize, and then read C 

the operational data for the triplet 
9930X/9930Y from the card side re- 
corded in the ODCP example. 

The indicated time delay is 28800 ys 
from 99 3 OX. 

Input the indicated time delay, then 28800 A 

press A. 

The indicated time delay is 49400 ys 
from 9930Y. 

Input the indicated time delay, then 49400 B 

press B. 

Press D to obtain the 1st solution D 

latitude . 

Press R/S to obtain the 1st solution R/S 

longitude . 

The first solution is: 42°44'57"N Latitude 

41 °07 ' 32"W Longitude 

Press E to obtain the 2nd solution e 

latitude . 

Press R/S to obtain the 2nd solution R/S 

longitude . 

The 2nd solution is: 27°00’07”S Latitude 

102°27’12”E Longitude 



Output 

Data/Units 



2B800 

49400 

42.4457 

-41.0732 

- 27.0007 

102.2712 
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3 , Program Listings 



Master Data Card Program 



0 ili ID 


76 LBL 


IJ 5 IJ 


— * c- 

1 J 


00 1 


4 2 3TG 


Ij 5 1 


Ij 2 2 


ft Tj V 


72 ST-^ 


052 


1 J 9 9 


f! !'■{ : 


Cl IJ 0 0 


053 


Ij 8 Cl 


004 


97 ns 2 


05 R 


93 


f\ iD ^ 


r; fi 0 ij 


Cl 5 5 


02 2 


iD L*i 0 


42 STO 


058 


Cl F. 8 


fi i"i 7 


0 i'l 0 


0 5 7 


35 


1 J 0 F: 


22 IH'v' 


1 J 5 1 J 


4 2 S 1 D 


1 J 0 9* 


58 FI 7! 


059 


2 1 J 2 1 J 


ij_L0 


04 4 


1 J 8 0 


C cr — 


01 1 


b b i“’ n L» 


08 1 


22 INV 


0 1 2 


b F. PR L* 


Ci 1 2 


30 1 RN 


0 1 C:: 


ft b p n 1 J 


Cl 8 C-i 


i' 2 b ! 


0 1 4 


Rb PRU 


1 J F. 4 


0 IJ 0 1 J 


0 1 5 


66 PflU 


Cl F. 5 


Ci iJ iJ 


Ij 1 b 


9 b hi h! i 


Cl 8 8 


92 RTN 


0 1 7 


92 Pi N 


i_i F. 7 


78 LBL 


Li 1 8 


76 LBL 


0 8 Cl 


12 B 


0 1 9 


1 ,■■■ B ' 


0 8 9 


88 II MS 


0 2 Ct 


42 STG 


0 7 Cl 


42 STO 


02 1 


2 9 2 9 


07 1 




022 


92 RTH 


0 7 2 


ijQ 0 


023 


1' 6 L D L. 


iJ 7 Cl 


92 R i N 


ij 2 H- 


16 PJ 


074 


76 LBL 


025 


42 SlO 


i-i *r« cr 

IJ >• •_> 


14 ii 


ij 2 R 


ill C' 2 c* 


0 7 F: 


88 DM 8 


0 2 7 


92 RTH 


077 


42 STG 


028 


7 6 L B L 


ij 7 Cl 


24 24 


029 


11 H 


079 


00 0 


0 0 




IJ 8 0 


92 K i H 


0 1 


iD 2 2 


08 i 


{'' 6’ LBL 


0 3 2 


ijS 3 


1 J 8 2 


1 5 E 


033 


6' 1 i-j i D 


0 Cl 3 


:3 RCL 


0 34 


85 + 


084 


2 jl 2 2 


Cl Ci 5 


76 LBL 


0 85 


•I« cr 


0 3 b 


1 3 C 


0 8 8 


43 RCL 


- - — . 


l, .1 ■. 'T- 


C' 


p; 2 5 


IJ i« r' 


-Z- iZZ : *’i ■* ; 


L* C' ! 




0 3 £' 


1 J 2 2 


088 


95 


iJ 3 9 


fi 5 5 


ij Cl 9 


55 ^ 


0 4 ij 


78 LBL 


0 9 IJ 


02 2 


04 1 


85 


09 1 


O cr _ 


042 


4 2 •!= i D 


ij 9 2 


42 STD 


-1 A “i 

IJ I- 


i J ij Ij Ij 


f ! Q ;3 


0 1 j Ci 0 


044 


2 ^ 1 


094 


4 3 RCL 


iJ 4 5 


8 IJ ji i 1 


ij 9 5 


25 2 5 


IJ 4 8 


3u ; RN 


Ij 9 F: 


i L- "" 


04 7 


J cr 


0 9 7 


4 i;l K L- L 


} J 4 !r' 


cr. “i 


Cl 9 8 


2 Cl 2 3 


1 J 4 9 


■i i 1 


ij 9 9 


‘^5 - 



1 00 
101 
1 02 
102 
1 04 
105 

I Ob 
107 
lOS 
1 09 

I I 0 
1 1 1 
i i Ll 

1 1 3 
114 

I 1 5 

i 16 

I I 

ns 

119 

120 
121 
I J 

123 
1 2 4 
i 25 

1 27 
129 

129 

130 

131 

1 ■■ ■> 
J. w 

133 

134 

135 
1 3 6 

137 

138 

139 

140 

141 

142 

143 

144 

145 
1 46 
1 47 

148 

149 



cr cr 



02 

95 

42 
01 

24 

75 

43 



c- T n 

Oi 

PCL 

■'-■a 



RC 



i 

i— 



9 n 
42 
02 
55 

02 
Q ^ 

42 
0 : 




t' j 

4 -2« 




42 

04 







1 ■ 



n 11 



G 



STD 



O T ki 
■-'ill 



ROL 

CDS 



R C L 

CDS 

R i_- L 
0 1 
C; T \l 



Sin 



04 






i 



STD 



F! L- L 

I’i 0 



15 



Master Data Card Program (continued) 



150 


4 3 


POL 


2 0 1 J 


4 3 


RCL 


250 


65 




151 


01 


fj 1 


201 


13 


IS 


251 


cr -J 




150 




CTJ6 


2 Ij 2 


3 3: 


S I M 


•2' 


Cl4 


•■t 


15 3 


95 


“ 


202 


65 




cr 

i2. '2= 


” cr 


- 


154 






204 




■n’ 


254 


43 


PCL 


155 


65 




2 Cl 3 


cr cr 


- 


255 


n 6 


fi 


156 


02 




206 


lii 


1 


256 


54 




15? 


cr cr 
■J J 


- 


207 


0 S’ 


Si 


257 


95 


= 


156 


cr ■“ 




2 ij S' 


00 


ij 


cr 

•2. J C 


65 




156 


0 1 


1 


2 IJ 9 


95 


= 


259 


43 


RCL 


160 


75 


- 


210 


42 


S- 0 


2 6 0 


20 


20 


161 


4 3 


POL 


211 


05 


03 


261 


65 




162 


04 


04 


212 


65 




262 


4 3 


RCL 


163 


95 




213 


43 


RCL 


2 S' -2' 


05 


05 


164 


42 


6TD 


214 


06 


Ij Si 


264 


65 




165 


0 b 


0 b 


215 


75 


- 


265 


43 


RCL 


166 


43 


POL 


2 1 6 


43 


RCL 


2 6 6 


02 


02 


lb? 


01 


01 


217 


09 


0 9 


2t’7 


3 0 


TfiH 


166 


t*: 


:3 1 H 


2 1 6 


95 


= 


2 6 3 


65 




169 


65 




219 


65 




269 


04 


4 


170 


43 


RCL 


2 2 ij 


43 


RCL 


270 


05 


5 


171 


0 0 


0 Cl 


221 


20 


20 


271 


55 


- 


172 


3 9 


CDS 


2 2 2 


55 




4 1 2 


69 


it 


176 


95 


= 


2 2 


04 


4 


*2. i‘ -2' 


94 


-i- / - 


174 


33 


i- •! ^ 


224 


95 




274 


65 


T 


175 


65 




cr 


42 


STD 


275 


43 


RCL 


1 , 6 


02 




226 


1 1 


1 1 


276 


02 


02 


177 


55 


- 


•~1 J “ 


94 


-r .•••■ - 


•2. !‘ i' 


95 


=: 


178 


43 


RCL 


2 2 S' 


.“r cr 
J 


-r 


i2. i' 0 


c c 

.J _l 




179 


04 


04 


2 9 


43 


RCL 


2 79 


02 


V 


160 


95 


= 


2 K 0 


n? 


n-. 


2S 0 


95 


= 


161 


42 


STD 


2-1 


95 


= 


281 


42 


5 TO 


162 


0 r" 


0? 


O 


65 






12 


12 


163 


i"; cr 

l-l .J 


i- 


-I' 


4 3 


RCL 




Z‘ 'Z- 


5 I H 


164 


43 


RCL 


234 


13 


13 


264 


65 




165 


06 


06 


“t cr 
«il -2' • J 


•2' C' 


6 I H 


265 


43 


RCL 


166 


95 


=: 


2 Si 6 


95 


= 


266 


fi M 


fi fi 


167 


42 


STD 


i2- -2* !' 


42 


STD 


•2. r 


:3 9 


CD6 


186 


1 IK 


i IK 


li 2> C' 


26 


2 6 


'2' 0 '2* 


95 


= 


169 


43 


RCL 


2 "i 


43 


RCL 


289 


i4- 


‘ r T 


1 9U 


Ub 


Ub 


240 


09 


09 


2 9 1 j 


43 


RCL 


191 


— ir 


- 


241 


—1 ^ 

.J 


- 


291 


12 




1 92 


4 2 


RCL 


242 


cr -i 

J J: 




2 K 2 


ij s 


C06 


1 5^3 


0^‘ 


ij 


243 


0 1 


1 


2 9 ij 


63 




1 7-t 


r. cr 

r 


- 


244 




- 


294 


4 3 


RCL 


1 95 


4 2 


S”u 


245 


' : K 




293 


0 1 


01 


196 


U r 


ij 


2 


63 




2 9 S 


3- . 3 ; 


6Irl 


1 97 


•: “ 


rr= 1 


24"^ 


•i - 3 . 


RC L 


297 


‘1: ~ 




196 


1 


1 :j 


24 S 


04 


04 


2 S: 


’ j 


i- .•••■ - 


199 


cr cr 


- 


24 9 


54 




2 9 9 


~= 7' 


INV 



16 



Master Data Card Program (continued) 
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Operational Data Program Card (continued) 
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Position Determination Program (continued) 
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Position Determination Program (continued) 
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D. 



Loran-C Fixing Algorithms 



The development of the Loran fixing algorithms in this 
report is presented in more detail in a companion report 
[Ref. 3] and will not be repeated here. 

The basic Loran-C equation [Ref. 4] can be written as 

T = [Tg + p(Tg)l - [T^ + + iTg + P(Tj^)] + ^ d) 

where 

T is the "indicated time difference" in microseconds, 
Tj^,Tg is the distance, in microseconds, from the master 
and the slave to the receiver, respectively, 

Tq is the distance, in microseconds, between the 

D 

master and the slave, 

6 is the assigned coding delay, in microseconds, and 

p(T) is the secondary phase correction, in microseconds, 
for an all sea water path of length T. 

The quantity 

At = [Tg + p(Tg)] + 6 

is a constant for each master/slave pair. The following 
World Geodetic System 1972 (WGS 72) values have been adopted 
for Loran-C navigation [Ref. 4]: 
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Vq - 299792458 meters/second is the velocity of light 
in free space, 

ii = 1.000 !38 is the index of refraction of the surface 
ot th-i earth for standard atmosphere and lOOkHz 
elect “omagentic waves, 

a^ = 6378135.00 meters is the equatorial radius of the 
earth, and 

f = 1/298.26 is the flattening factor (1 - b/a^, where 
b is the polar radius) of the earth. 

Accurate formulas for computing the secondary phase 
correction p(T) are contained in Reference 4, but for use 
with the handheld calculator the following linear approximation 
[Ref. 3] will be used: 



where 

and 



p(T) = a^^ + ^2^ ' 



a^ = -0.321, 

= 0.000635. 



Using this approximation, it is possible to solve Equation 1 
for the quantity Tg - We find 



Tg - = (T - At)/(1 + a^) 



( 2 ) 



On the surface of a sphere a hyperbolic line of position 
can be represented by the equation [Ref. 3, page 175] 
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tan r 



cos 2a - cos 2c 

sin 2c cos w + c sin 2a 



(3) 



where the origin of the coordinate system is at the prime focus 
of the spherical hyperbola, 2c is the spherical arc joining 
the foci, 2a is a constant for any one LOP, and r and uj 
are the spherical coordinates of a point on the LOP. If the 
base line of the coordinate system is the arc joining the foci 
then u) is the spherical polar angle from the base line to a 
point P on the LOP and r is the spherical polar distance 
(or arc) from the prime focus to P. Using the Loran system 
we take f, = +1 if the prime focus is at a master station 
and we take C = -1 if the prime focus is at a slave station. 

If we take v = v^/n to be the velocity of lOOkHz electro- 
magnetic radiation of the earth's surface then 



or , using Eq . ( 2) , 



2a = v(Tg - T„)/a^ , 



2a = (T - At)/a , 

It 



(4) 



where 



a (1 + a ) 

a = — 7 — — = 21295.87 us . 

P v_/n 



The baseline between master and slave can be obtained from 



2c = V Tg/a^ . (5) 

Here 2c is computed by program card 1 (preparation of master 
data cards) using the algorithm in Section E. 
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Consicior u Loran-C triplet with master stations colocated. 
Let ^ and ^ denote the azimuth angles of slave 1 (S^) 

and slave 2 J^^-spectively , measured from North toward the 

East from the master stations (M) (see Fig. 3) . Further, let 
u and r denote the azimuth and spherical polar arc (distance) 
of the receiver (R) from M. For this geometry, Eq. (3) can 
be written as 



where 



tan r . 
1 



B . 
1 



C . cos (a~C • ) + A . 

X X 



( 6 ) 





> 

II 


c, . sin 2a . 
1 1 








B . = 

X 


cos 2a. ~ cos 
1 


2c . 

X 




and 


C. = 

X 


sin 2c . 

1 






for 


. . . th 

the X — 


Loran pair, i = 


1,2. 


Since r = r, = tan r 



can be eliminated in Eq . (6) . The resulting equation can be 

rewritten as 



C cos u+Ssina=K, (7) 



where 





C = B^C^ 


cos 


- 


COS 






S = B^C2 


sin f, 2 


- 


sin 




and 


K = B^A^ 


®1^2 


• 
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N 




Figure 3. Geometry of a Loran Triplet and a Receiver. 
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It wu ciut ine i 


> 0 and , by the equations 


and 


. cos 1 - C, 

(8) 

• sin , = S, 


then 


P - /c^ + , 


and 


y = qatn(S,C) 



Here the function (.jatn(y,x) is ttie arctangent of y/x 
adjusted for the proper quadrant according to the signs of x 
and y. A compact form of this function is 



qatn(y,x) = tan ^ 


^ + 11 1 { x < 0 ? ) 

X + 10 t(x - 0?) 



where 



t(z) = 1 when 


z is true 



and 



t ( z) =0 when 


z is false. 


When convenient we will 


use the notation qatn(y/x) inter- 



changeably with qatn(y,x) . The qatn function is equivalent 
to the polar angle obtained using the rectangular to polar 
conversion function on the TI-59. 

Now substitute Eq . (8) into Eq . (7) and solve for 

a = Y + cos ^ ( K /p ) ( 9 ) 
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to obtain the azimuth angle a of the two points of intersection 
of the lop's. Finally we obtain a value for r by substituting 
each u into Eg. (5) . We find that 



r 



gain 



B . 

1 

C . cos (u-f, . ) + A . 
1 1 1 



for i = 1 or 2 . 



The distance and a 
converted into the 
positions of R. 

The fixing 
solution algorithm 



zimuth from M or the triplet vertex can be 
latitude and longitude of the two possible 

algorithm then uses u and r in the dived 
of spheroidal geodesy (Section F) . 
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E. 



'I'lie Reverse ( 1 nverse) Solution Algorithm 



This /•< solution algorithm is a modification of the 

first order in flc:ttening (f) algorithm given by Thomas (Ref. 5, 
pp. 8-10). Tliomas ' notation has been followed as closely as 
possible for ease of comparison of the algorithms. The gatn 
function is defined in Section D. West longitudes (\) and 
South latitudes ( >1 ) are negative. We are given the points 
Pj^ ( i|) / A ^ , P2(tl>2f^2^ spheroid and are to find the 

distance S between the points and the forward and back 
azimuths, '^12 ‘^21* Given quantities are A ^ 

and A^. No assumptions about the relative location of 
and P 2 are required. The modified vevei-^e solution algorithm 
is : 

= tan tail 4,^), i = 1, 2, 



(6, + e_) / 2 , Ab = 

12 m 


(e^-ej)/2. 


AA = 


^2 


- S' 


Lk/2, H = cos^A0 


. 2 

- sin 0 = 


2 

cos 


0 


. 2 

- sin AB = cos 0, cos 0 



m m m m 12 



2 2 2 2 
L = sin + H sin AA^ = sin (d/2) , 1 - L = cos (d/2) , 

- 1 2 2 
d = cos (1 - 2L), U = 2 sin 6 cos A9 /(I - L) , 

m m 

2 2 

V = 2 sin A6 cos 0 /L, X=U+V, Y=U-V, 
mm ' 

T = d/sin d, 6 ^d = f(TX - Y)/4, S = a^(T-6jd) sin d, 

F = 2{y - (1 - 2L) (4 - X) ] , G = fT/2, 

Q = -(FG tan AA)/4, AA ' = (AA + Q) /2 

m ^ 
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qcitn(-iiin AO cos AA * , cos 0 sin AA ' ) , 
m m m m 

qutn(cos AO cos AA * , sin sin AA * ) , 

in m rn m 

S " q- -21 = S - q • 

This reverse solution algorithm is used by program card 1 

(preparation of master data cards) to compute the baseline 

distance 2c and the azimuths and between the 

MS SM 

master and slave stations of a Loran pair. 

Details of the modifications made to Thomas' algorithm 
are contained in Reference 3. 
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F. 



The Direct Solution Algorithm 



This /■/'(•/ solution alyorithm is a modification of 
the first order an flattening (f) algorithm given by Thomas 
(Ref. 5, pp. 7-8] . Thomas' notation has been followed as 
closely as possible for ease of comparison of the algorithms. 
The qatn function is defiraed in Section D. West loiagitudes 
and South latitudes are negative. We are given the point 
^l^'^l'^l^ on the spherioid, where ^ are the geodetic 

latitude and longtiude (geographic coordinates); the forward 
azimuth u ^2 distance S to a secorad point 

and from these we are to find the geographic coordinates 

A 2 the back azimuth '^21' given quantities are 

4) , Aw a. „ and S. No assumptions about the relative loca- 

X i. X 

tion of and P 2 are required. The modified dlvect 

solution algorithm is: 

= tan ^[(1-f) tan M = cos 0^^ sin 

2 

N = cos cos ^^ 2 * = f ( i - M ) /4 , 

D = 1 - 2c^ c^M, P = = qatn(N^ sin 0^) 

d = S/(a^D), u = 2(0^--d), W = 1 - 2P cos u, 

V = cos(u + d) , y = 2PVW sin d, Ao = d-Y , 

“ qatn[-M, - (N cos Ao - sin sin Ao ) ] , 

K = (l“f) [m"^ + (N cos Ao - sin sin Ao)^]^"^^ , 

't* ^ = tan ^ [ (sin 6 ^ cos Ao x N sin Ao ) /K ) , 



An 


= qatn (sin Ao 


s i n 


‘^12' 


cos fos Ao - sin 6^ sin Ao cos 


12 ' ' 


li 


II 

0 

> 

> 

H 


Ari 


^ H, 


A-, = >., + AA. 
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This direct solution algorithm is used by program 
card 3 ( improvt;d fix program) to compute the latitude and 
longitude of the receiver using the azimuth and range of the 
receiver from the Loran triplet vertex. 

Details of the modifications made to Thomas' algoritlim 
are contained in Reference 3. 
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G. Some Results 



The following procedure was used to test the fixing 
algorithm: A more accurate version of the reverse solution 

algorithm (see Reference 3) was used to compute time delays 
at selected positions from Loran-C pairs. These computed 
time delays do not depend upon a fixing algorithm. The 
fixing algorithm was tested by determining the distance 
betv/een the fix and the positions used to compute the time 
delays. Typical results are presented in Tables I through 
IV. As can be seen, all differences (A) are within 1 n.mi. 
The greatest difference (A = 0.92 n.mi. in Table IV) was 
obtained at a position located 2270 n.mi. from the triplet 
vertex. 
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Table I. 9940 X/ 9940 Y 



Position 
Lat Long 


Calculated 

9940X 


Time Delay 
9940Y 


TI-59 
Lat (N) 


Fix 

Long (W) 


A 

n.mi 


24°N 


122®W 


27726.19 


40912.76 


23^59' 


155.. 


122^00' 


01” 


0.08 


26 


122 


27715.97 


40998.39 


25°59 


.57.. 


122^00' 


'01” 


0.05 


28 


122 


27702.41 


41117.84 


28°00 


'00” 


122^00’ 


'00” 


0.00 


30 


122 


27683.53 


41291.85 


29059 


'59” 


122^00' 


■00" 


0.02 


32 


122 


27655.47 


41555.46 


32°00 


*00” 


122^00’ 


'00” 


0.00 


34 


122 


27609.63 


41959.57 


U) 

6 

0 


'00” 


122^00' 


'00” 


0.00 


36 


122 


27523.56 


42544.11 


36^00 


*00” 


121^59' 


*59” 


0.01 


38 


122 


27334.61 


43248.22 


0 

0 

0 

CO 

m 


*00” 


122°59' 


*59” 


0.01 



Table II. 9940 Y/ 9940 W 



Position 


Calculated 


Time Delay 


TI-59 


Fix 


A 


Lat 


Long 


9940Y 


9940W 


Lat (N) 


Long (W) 


n.mi 


37 


122°W 


42892.86 


16257.23 


36°59'59" 


122°01'01" 


0.02 


37 


125 


43056.68 


15765.13 


U) 

0 

0 

0 

0 

0 


125"00'00" 


0.00 


37 


128 


43137.78 


15327.12 


37°00'00" 


128°00'00" 


0.00 


37 


131 


43191.10 


14970.77 


U) 

0 

0 

0 

0 

0 


131°00'00" 


0.00 


37 


134 


43232.38 


14683.74 


b 

0 

0 

0 

0 

CO 


134®00'00" 


0.00 


37 


137 


43267.42 


14449.40 


37°00'00" 


137°00'00" 


0.00 


37 


140 


43298.80 


14254.02 


U) 

0 

0 

0 

0 

0 


140°00'00" 


0.01 


37 


143 


43327.85 


14087.43 


37°00'01" 


142'’59'59" 


0.02 
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Table III. 7930Z/9930X* 



Position 
Lat Long 


Calculated 

7930Z 


Time Delay 
9930X 


TI-59 
Lat (N) 


Fix 

Long (W) 


A 

n .mi 




30°W 


52437.86 


28451.72 


60°00'03” 


29059 . 32 »' 


0.24 


58 


35 


51960.93 


28391.50 


b 

o 

o 

o 

o 

00 

m 


34°5?*48" 


0.09 


56 


40 


50992.37 


28359.15 


55059159 .. 


39059 . 54 .. 


0.06 


54 


45 


49292.46 


28370.85 


53°59*59" 


44059 . 57 .. 


0.03 


52 


50 


47165.60 


28490.64 


52°00*00" 


49°59'57’‘ 


0.01 


50 


55 


45236.59 


29070.48 


50°00*00" 


55°00'00" 


0.00 


48 


60 


44505.60 


30991.94 


0 

0 

0 

0 

0 

00 


0 

0 

0 

0 

0 

0 


0.00 


46 


65 


44475.70 


33697.14 


46°00*00” 


65°00*00” 


0.00 


44 


70 


44588.91 


36567.42 


43059 . 59 .. 


69“59'59” 

‘ 


0.02 



Table IV. 9930W/9930X* 



Position 
Lat Long 


Calculated 

9930W 


Time Delay 
9930X 


TI -59 ] 

Lat (N) 


Fix 

Long (W) 


A 

n .mi 


9^N 


47^W 


13058.04 


36466.46 


3059159 .. 


46°59’22” 


0.92 


12 


52 


12984.71 


37288.35 


11^59*34” 


51 °59’ 37^' 


0.57 


15 


57 


12898.73 


38267.58 


14059 . 45 .. 


56°59*47” 


0.28 


18 


62 


12793.91 


39431.32 


17059 . 52*» 


61°59'54” 


0.16 


21 


67 


12656.52 


40794.36 


20°59*56'‘ 


66°59’57'' 


0.08 


24 


72 


12451.30 


42330.55 


23059 . 59 .. 


71059 . 59 .. 


0.02 


27 


77 


12097.12 


43876.62 


t — 1 

0 

0 

0 

0 

CN 


77°00'00” 


0.02 


30 


82 


12973.95 


44768.53 


30°00'01" 


82°00'06” 


D.09 



*The 9930 system was disestablished 30 September 1979. It is included here 
for illustration only. 
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APPENDIX. Loran-C Station Parameters 



The following list contains the pertinent parameters 
for each Loran-C station pair. This list was compiled from 
data in Reference 4. Each row contains the following inform- 
ation : 

1. The Loran-C station pair designator 

2. At, the sum of the coding delay plus one way base- 
line time, including the secondary phase correction 
for an all seawater path, in micrseconds. 

3. The master station latitude. 

4. The master station longitude. 

5. The slave station latitude. 

6. The slave station longitude. 

In this report, the convention that a minus (-) sign denotes 
South latitudes and West longitudes has been followed 
throughout Some users may wish to reverse the longitude 
convention so that (-) longitudes are East rather than West. 
This may be accomplished by consistently changing the sign 
of the longitude everywhere that it is required for data 
input; no program modifications are required. If the sign 
convention for longitude is reversed, care should be taken 
that the signs of all longitudes in the list are reversed. 

In columns 3 through 6 the latitudes and longitudes appear 
to be in decimal form, but the actual format is DDD.MMSSFF 
where 
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DDD designates degrees, 

MM designates minutes, 

SS designates seconds, and 

FF designates hundredths of seconds. 
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Coveracje of Loran-C Systems 



Station 


Location 


4990 


Central Pacific 


5930 


East Coast, Canada 


5990 


West Coast, Canada 


7930 


North Atlantic 


7960 


Gulf of Alaska 


7970 


Norwegian Sea 


7980 


Southeast U.S.A. 


7990 


Mediterranean Sea 


8970 


Great Lakes 


9930 


East Coast, U.S.A. 


9940 


West Coast, U.S.A. 


9960 


Northeast U.S.A. 


9970 


Northwest Pacific 


9990 


North Pacific 
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